gsem  (dichotcont <-  policyrepdrift i.urban age total_child i.maritcat ib3.educat   i.wealth_cat  ) (pt1_tri1_estim_induced <- policyrepdrift i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol  dichotcont famplan_dis11 )if age >19 & age <41, poisson

estat ic


nlcom _b[pt3:modern]*_b[modern:policyrepdrift]
nlcom (_b[pt3:modern]*_b[modern:policyrepdrift]+_b[pt3:policyrepdrift])
nlcom _b[pt3:modern]*_b[modern:policyrepdrift]/ (_b[pt3:modern]*_b[modern:policyrepdrift]+_b[pt3:policyrepdrift])




gsem  (planvis <-  policyrepdrift i.urban age total_child i.maritcat ib3.educat   i.wealth_cat  ) (pt3 <- planvis policyrepdrift i.urban age total_child i.maritcat ib3.educat  i.wealth_cat )if age >19 & age <41, poisson

estat ic


nlcom _b[pt3:planvis]*_b[planvis:policyrepdrift]
nlcom (_b[pt3:planvis]*_b[planvis:policyrepdrift]+_b[pt3:policyrepdrift])
nlcom _b[pt3:planvis]*_b[planvis:policyrepdrift]/ (_b[pt3:planvis]*_b[planvis:policyrepdrift]+_b[pt3:policyrepdrift])



melogit  pt1_tri1_contr     law i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic
melogit  pt1_tri1_contr     policyrep i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic
melogit  pt1_tri1_contr         policyrepdrift i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic

melogit  pt1_tri1_unwanted     law i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic
melogit  pt1_tri1_unwanted     policyrep i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic
melogit  pt1_tri1_unwanted         policyrepdrift i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic


melogit  pt1_tri1_estim_induced     law i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic
melogit  pt1_tri1_estim_induced     policyrep i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic
melogit  pt1_tri1_estim_induced         policyrepdrift i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol ib2.REGION i.current_type_cat famplan_dis11 if age >19 & age <41  ||int_year: || cluster_num:, or
estat ic




****Contraceptive resuls for SEM***
gsem  (dichotcont <- policyrepdrift famplan_dis11 law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year ) (pt1_tri1_estim_induced <- policyrepdrift  i.edu_level i.age_group ib4.maritst_ca dichotcont famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit) 
estat ic

nlcom _b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrepdrift]
nlcom (_b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrepdrift]+_b[pt1_tri1_estim_induced:policyrepdrift])
nlcom _b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrepdrift]/ (_b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrepdrift]+_b[pt1_tri1_estim_induced:policyrepdrift])

gsem  (dichotcont <- policyrep famplan_dis11 law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year ) (pt1_tri1_estim_induced <- policyrep  i.edu_level i.age_group ib4.maritst_ca dichotcont famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit) 
estat ic

nlcom _b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrep]
nlcom (_b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrep]+_b[pt1_tri1_estim_induced:policyrep])
nlcom _b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrep]/ (_b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:policyrep]+_b[pt1_tri1_estim_induced:policyrep])

gsem  (dichotcont <- law famplan_dis11 law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year) (pt1_tri1_estim_induced <- law  i.edu_level i.age_group ib4.maritst_ca dichotcont famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit)

nlcom _b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:law]
nlcom (_b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:law]+_b[pt1_tri1_estim_induced:law])
nlcom _b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:law]/ (_b[pt1_tri1_estim_induced:dichotcont]*_b[dichotcont:law]+_b[pt1_tri1_estim_induced:law])

****Family conversation resuls for SEM***


gsem  (famplan_dis11 <- law  dichotcont i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year) (pt1_tri1_estim_induced <- law  i.edu_level i.age_group ib4.maritst_ca dichotcont famplan_dis11 i.urban total_child domviol i.regionsem i.int_year  )if age >19 & age <41,  family(bernoulli) link(probit)

nlcom _b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:law]
nlcom (_b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:law]+_b[pt1_tri1_estim_induced:law])
nlcom _b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:law]/ (_b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:law]+_b[pt1_tri1_estim_induced:law])

gsem  (famplan_dis11 <- policyrep dichotcont law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year ) (pt1_tri1_estim_induced <- policyrep  i.edu_level i.age_group ib4.maritst_ca dichotcont famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit) 

nlcom _b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrep]
nlcom (_b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrep]+_b[pt1_tri1_estim_induced:policyrep])
nlcom _b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrep]/ (_b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrep]+_b[pt1_tri1_estim_induced:policyrep])

gsem  (famplan_dis11 <- policyrepdrift dichotcont law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year) (pt1_tri1_estim_induced <- policyrepdrift  i.edu_level i.age_group ib4.maritst_ca dichotcont famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit) 

nlcom _b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrepdrift]
nlcom (_b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrepdrift]+_b[pt1_tri1_estim_induced:policyrepdrift])
nlcom _b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrepdrift]/ (_b[pt1_tri1_estim_induced:famplan_dis11]*_b[famplan_dis11:policyrepdrift]+_b[pt1_tri1_estim_induced:policyrepdrift])



***Access**
gsem  (dichorigcontrac <- law  dichotcont i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year) (pt1_tri1_estim_induced <- law  i.edu_level i.age_group ib4.maritst_ca dichotcont dichorigcontrac famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit)

nlcom _b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:law]
nlcom (_b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:law]+_b[pt1_tri1_estim_induced:law])
nlcom _b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:law]/ (_b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:law]+_b[pt1_tri1_estim_induced:law])

gsem  (dichorigcontrac <- policyrep dichotcont law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year) (pt1_tri1_estim_induced <- policyrep  i.edu_level i.age_group ib4.maritst_ca dichorigcontrac famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit) 

nlcom _b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrep]
nlcom (_b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrep]+_b[pt1_tri1_estim_induced:policyrep])
nlcom _b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrep]/ (_b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrep]+_b[pt1_tri1_estim_induced:policyrep])

gsem  (dichorigcontrac <- policyrepdrift dichotcont law  i.edu_level i.age_group ib4.maritst_cat i.urban total_child i.domviol i.regionsem i.int_year) (pt1_tri1_estim_induced <- policyrepdrift  i.edu_level i.age_group ib4.maritst_ca dichorigcontrac famplan_dis11 i.urban total_child domviol i.regionsem i.int_year)if age >19 & age <41,  family(bernoulli) link(probit) 

nlcom _b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrepdrift]
nlcom (_b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrepdrift]+_b[pt1_tri1_estim_induced:policyrepdrift])
nlcom _b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrepdrift]/ (_b[pt1_tri1_estim_induced:dichorigcontrac]*_b[dichorigcontrac:policyrepdrift]+_b[pt1_tri1_estim_induced:policyrepdrift])




sem  (famplan_dis11 <- policyrepdrift  ) (pt1_tri1_estim_induced <- policyrepdrift    )if age >19 & age <41,   
sem  (famplan_dis11 <- policyrep  ) (pt1_tri1_estim_induced <- policyrep    )if age >19 & age <41,   
sem  (famplan_dis11 <- law  ) (pt1_tri1_estim_induced <- law    )if age >19 & age <41,   


sem  (dichorigcontrac <- policyrepdrift  ) (pt1_tri1_estim_induced <- policyrepdrift    )if age >19 & age <41,   
sem  (dichorigcontrac <- policyrep  ) (pt1_tri1_estim_induced <- policyrep    )if age >19 & age <41,   
sem  (dichorigcontrac <- law  ) (pt1_tri1_estim_induced <- law    )if age >19 & age <41,   


sem  (dichotcont <- policyrepdrift  ) (pt1_tri1_estim_induced <- policyrepdrift    )if age >19 & age <41,   
sem  (dichotcont <- policyrep  ) (pt1_tri1_estim_induced <- policyrep    )if age >19 & age <41,   
sem  (dichotcont <- law  ) (pt1_tri1_estim_induced <- law    )if age >19 & age <41,   
